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Abstract. In image processing, it can be a useful pre-processing step to smooth away small 
structures, such as noise or unimportant details, while retaining the overall structure of the image 
by keeping edges, which separate objects, sharp. Typically this edge-preserving smoothing process 
is achieved using edge-aware filters. However such filters may preserve unwanted small structures as 
well if they contain edges. In this work we present a novel framework for edge-preserving smoothing 
which separates the process into two different steps: First the image is smoothed using a blurring 
filter and in the second step the important edges are restored using a guided edge-aware filter. The 
presented method proves to deliver very good results, compared to state-of-the-art edge-preserving 
smoothing filters, especially at removing unwanted small structures. Furthermore it is very versatile 
and can easily be adapted to different fields of applications while at the same time being very fast 
to compute and therefore well-suited for real time applications. 


1 Introduction 


Edge-preserving smoothing is an approach to smooth across boundaries between regions of similar ap¬ 
pearance while retaining strong edges. It is a technique to remove noise, weak edges and small details 
whereas the overall structure of the image should not be lost. Recent edge-preserving smoothing tech¬ 
niques, such as the joint bilateral filter [7], the guided filter [ 3 ] and the rolling guidance filter 10 , use a 
so called guidance image to perform this operation. 

In this work we propose a novel framework for edge-preserving smoothing. Like the rolling guidance 
filter proposed in 10 it aims not only for edge-aware filtering, but also for scale-aware filtering. Our 


model splits the tasks of smoothing and edge preservation by using two different filters: a smoothing 
filter and a guided edge-aware filter. In the first step the smoothing filter removes small structures from 
the input image, but it also blurs strong edges. The second step tries to restore those smoothed strong 
edges iteratively by applying the guided edge-aware filter. 

In Section [2] of this work we present a short review of recent edge-preserving filters which are most 
related to our work. Section [3] describes our edge-preserving smoothing model in detail and gives some 
possible choices for the smoothing and the edge-aware filters. It further introduces a simple edge-aware 
filter, by using an approach similar to [8] and discusses its properties. The results of our experiments 
and possible applications of this new framework are presented in Section [4] and finally Section [5] gives a 
short summary and presents possible prospective work on the topics discussed. 


2 Related work 

In the past, edge-aware filters were a popular approach to edge-preserving smoothing. These are filters 
that detect large intensity differences between neighbouring pixels and classify such differences as edges 
to be preserved (strong edges). Low differences (weak edges ) are treated as unimportant structures 
and smoothed away. Early approaches include for example anisotropic diffusion 6 or the symmetric 
nearest neighbour filter ( 2 ], an example for a more recent approach is the guided filter ( 3 ]. Such filters 
normally ignore the size of the structures in an image, as small structures can contain strong edges as 
well. However, such small structures are often unwanted details that we want to smooth away. A filter 
that takes into account the size of an edge is called a scale-aware filter. Edge-aware filters are therefore 
usually not scale-aware. 

Whereas edge-aware filters only detect edges with large intensity difference, the bilateral filter 9 , being 
a combination of a Gaussian in the spatial domain and a Gaussian in the range domain, smoothes away 
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small structures and therefore detects only strong edges larger than a certain scale, depending on its 
parameters. While the bilateral filter generates satisfying output images, it is very slow to compute. 
Therefore different optimizations have been proposed to speed up the computation (e.g. , but 

these are only approximations to the original bilateral filter, suffering from a loss in quality or needing 
highly optimized code running on GPU’s, especially if aiming for real-time performance. To further 
enhance the bilateral filter, the joint or cross bilateral filter was introduced [ 7 ], which uses a guidance 
image to better remove noise or unwanted structures in an image. 

Using this joint bilateral filter and exploiting its scale-awareness, the rolling guidance filter introduced 
in 10 iteratively computes a guidance image to apply the joint bilateral filter (or generally any edge¬ 
preserving smoothing filter) to the input image. It is the work that is most related to ours, as it aims 
for scale-aware filtering by blurring away small structures as well. 


3 Smooth and iteratively Restore 

While applying only edge-aware filters leads to good results concerning strong edges, they lack at remov¬ 
ing small structures from the image, because those small structures often contain strong edges themselves. 
Therefore it is necessary to remove those small structures before applying an edge-aware filter, so that 
the edge-aware filter can smooth the image while preserving strong edges. To achieve this, we first apply 
a different filter which removes small structures, but which doesn’t care much about strong edges, re¬ 
lieving the edge-aware filter from the task to remove small structures. This can be for example a spatial 
domain filter such as the Gaussian filter. However applying the edge-aware filter to the blurry output of 
the smoothing filter won’t result in any visually appealing images, it will most probably result in worse 
output than after applying the edge-aware filter on the original input. But we still have the original input 
image, which has the best edge information we can have. By extending the edge-aware filter to use the 
original input image as a guidance image, we obtain a guided edge-aware filter which has no problems 
filtering the blurry output of the smoothing filter. It chooses the strong edges from the original input 
image and accordingly restores the blurry output image of the smoothing filter. This whole process can 
be simply described as smooth and iteratively restore (SiR). 

3.1 The Algorithm 

In the following sections we present some possible choices for the smoothing and the guided edge-aware 
filter, but for now let us assume we have already chosen two such filters. Let us denote the smoothing 
filter with <S, the guided edge restoring filter 1Zi g using a guidance image Jg, our input image with / 
and our output image with O. Then our algorithm first performs a smoothing operation by applying S 
to the input image I. This smoothing operation removes small structures up to a certain scale. How 
much it removes depends on the parameters of S. Any smoothing filter that removes unwanted small 
edges and structures can be used here, e.g. the Gaussian filter or box filter. 

Now we have a smoothed image with already removed small structures, but the strong edges are smoothed 
as well. We therefore unsmooth those strong edges, while smoothed regions without strong edges are 
left as they are. We apply 1Zi g in an iterative manner using n iterations. This iterative restoring step 
cannot restore small edges: even though they are present in the guidance image, the necessary color 
information is not present any more in the blurred image. Therefore, when 7 Zi G tries to restore those 
small structures found in Jg, it has no color information available to do that and it will therefore not 
be possible to restore those edges. On the other hand, large edges are still visible in the blurred image, 
therefore 7 Zi G can still restore large edges present in Iq using the color information from the blurred 
image (see Figure [l] for a small example of how this works). 

Putting those two steps together, we get the following algorithm: 

Input: input image /, guidance image Jg, smoothing filter <S, guided edge-aware filter 7£/ G , 
number of iterations n 
Output: filtered image O 

Blur image / with smoothing filter S : O = S(I) 
for i = 1 ...n do 

1 filter O with guided edge-aware filter 7 Zi G : O = 7 Zj G (0) 

end 

Algorithm 1: Smooth and iteratively restore 
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It is easily implemented as soon as you have two filters which fulfill the requirements for S and 1Zi G and 
because of its very flexible nature, the algorithm can be easily adapted to the needs of your application. 




Figure 1: Example showing how the algorithm works. Left: original image (guidance image) with two strong 
edges, one is large (magnification on bottom), the other small (magnification on top). Center: blurred image, 
large structure still visible, small structure blurred away. Right: output of SiR after 5 iterations, large structure 
was restored successfully, while small structure could not be restored due to unsufficient color information in 
blurred image. 


3.2 Remove small structures 

To remove small structures from an input image / we apply a smoothing filter S. This can be for example 
a domain filter, such as a Gaussian filter or a box filter, but can generally be any filter that smoothes 
away small structures. For this first step, it is not important to choose an edge-aware filter, restoring 
strong edges is the goal of step 2, but you can of course choose your filter to be more complex (e.g. 
bilateral filter). We experimented in particular with two different smoothing filters: the Gaussian filter 
and the box filter, which can be applied several times to approximate a Gaussian filter. 


3.3 Restore strong edges 


After removing small structures we apply an edge-aware filter to restore strong edges; this is a filter 
that smooths differences in pixel intensities and not in pixel coordinates. For our framework, we need 
to choose the filter to be guided in that sense, that it uses a guidance image to perform the filtering 
operation, such as the guided filter [3] or the rolling guidance filter 10 . To restore the strong edges 
from the input image I we can just use I as guidance image, but for different applications (such as 
flash/no-flash denoising, see 0 ) we might want to use other guidance images. 

For the rest of this subsection we use following notation: the original input image (before applying S ) 
will be denoted /, the output image of the first step (this is <S(/)) will be denoted J and the output 
image of the whole algorithm will be denoted O. 

We now present three possible guided edge-aware filters to use in our framework. In the results section 
we show some output images generated with those filters and discuss some of their properties. 


3.3.1 2D Gaussian Range Filter A simple range filter, that is used for example in the bilateral 
filter, is the Gaussian function on the pixel intensities: 
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working in a neighbourhood N of (£, f/), that favours pixels (x,y) with similar color. It detects strong 
edges (this means (J(x,y) — J(x,y)) 2 is very high) and is therefore well suited to smooth while retaining 
strong edges. However to work on an already smoothed image J, we modify 0 to use a guidance image 
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with the normalization coefficient a = y)eN(x y) . Here we still filter the image 

J, but we use Iq as a guidance image. We therefore detect the strong edges of the guidance image Iq 
and restore them in the smoothed image J. 


3.3.2 Recognizing pixels in the same region The 2D Gaussian range filter has two drawbacks: 
first of all, it is very slow to compute (see Table [l] for a comparison). If our goal is a fast filter for 
real-time applications, this filter will not suit our needs, especially if we are going to use it iteratively. 
The second disadvantage is a bit more difficult to spot. What we actually want to achieve is some 
kind of segmentation: we want an edge-preserving algorithm to recognize different regions/objects in our 
image, to smooth intensity discrepancies inside those regions and to sharpen the borders between those 
regions. The Gaussian range filter detects regions only by comparing intensities of two pixels p and g, 
so it actually works with following definition of a region: two pixels p and q are in the same region when 
they have nearly the same intensity value. Here p and q are not necessarily direct neighbours. 

This is a somewhat simple definition of a region. Humans can distinguish different objects of similar 
color in an image as long as there is a visible border between those regions. The problem with the above 
definition of a region is, that two pixels p and q can be assigned to the same region even if they belong 
to different objects in the image, because this definition doesn’t take into account the possible occurence 
of a border between p and q. 

To overcome this issue we introduce paths. Given two pixels p and q in the image a pixel r is in between 
p and q if d(p, r) + d(r, q) < / • d(p, q) for a distance measure d(u,v) and a constant / > 1. Q We say 
a path from p to q is a finite sequence of not necessarily neighbouring pixels to, ..., tk for some k > 0, 
where to := p and tk := g, such that for every 0 < i < &, ti is in between t^_i and U+%. We call G, ..., 
tk— i the intermediate pixels between p and q. Furthermore we say that there is an edge between two 
consecutive pixels of the path ti and £* + i if there is an intensity difference (I(U) — I(U+ 1)) 1 2 > 0. This 
edge is called a weak edge if the difference is low and a strong edge if there is a large difference. Finally 
we say that there is a strong edge along the path from p to q if there exists an i > 0, such that there is 
a strong edge between ti and U+i* 

With these definitions we are now able to come up with a better definition for a region: two pixels p 
and q with a predefined path to, ..., tk lie in the same region if there is no strong edge along this path. 
Notice however that we didn’t specify how exactly this path should be chosen. Different paths from p 
to q may lead to different regions, so choosing a reasonable path is crucial for defining a useful region. 
Our first definition of a region is a special case of the second definition: for all pixels p and q it simply 
considers the trivial path consisting of only p and g. To see the difference compared to the general case 
with a path to, ..., tk, k > 1, we can distinguish four cases: 

I All ti have almost the same intensity, so both p and q he in the same region according to both 
definitions. 

II p and q have different values and there exists an i such that 0 < i < k where the intensity difference 
(. I(U ) — I{ti+ 1 )) 2 is very high. According to both definitions p and q he in different regions. 

III If the path from p to q is a monotonic ramp, i.e. if I{p) < I{t\) < ... < I{tk- 1 ) < /(g) or 
I(p) > I(ti) > ... > I(tk~ i) > /(g) such that there is a high intensity difference (I(p) — /(g )) 2 
between p and g but low intensity differences (/(£«) — I(U+ 1 )) 2 between consecutive points of the 
path, then the first definition would place p and g in different regions, whereas the second definition 
would recognize no strong edge but rather a smooth transition from p to q and would therefore place 
them in the same region. 

IV Finally if p and q have similar values but there exists an i such that 0 < i < k where (I(U) — I{ti+\)) 2 
is very high, then we have at least one strong edge along to, ..., tk- According to the first definition 
of a region, p and q belong to the same region but the second definition places them in two different 
regions. 


Only in the first two cases the first definition of a region leads to a desired classification. In the other 

as [TT] and [IV] as if it were [TJ 


two cases it doesn’t, because it treats III 


1 Notice that the choice d(u : v ) = \u x — v x 

q as between p and q 


| u y — v y | and f = 1 defines all pixels within the window spanned by p and 
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Using paths, we can better classify the relationship between p and q. The next filter uses a simple path 
with d(u,v ) = \u x — v x \ + \u y — v y \, f = 1 and one intermediate pixel t\. 


3.3.3 Separable Gaussian Range Filter We now proceed similarly to [8], who approximate the 2D 
bilateral filter by simply applying a ID bilateral filter first horizontally and then vertically. But instead 
of approximating the whole bilateral filter, we only approximate the Gaussian range filter ©• First we 
define two operators Ohj and O v j as follows: 

O hJ (J)(x,y) = a r 1 £ e-^ 1 ^- 1 ^ 2 J(x,y) (3) 

xEN 1 (x) 

O vJ (J)(x,y) = ap £ J(x,y) (4) 

y£N 2 (y) 

with normalization coefficients aq and and neighbourhoods N\ of x and N 2 of y. Those two oper¬ 
ators compute a Gaussian over the pixel intensities, but they do it only in one dimension: Ohj filters 
horizontally and 0 Vl vertically. The combination Oh,i(O v j(J)) of these two operators leads to a filter 
that can be written in a closed form as 

0(x,y) = a - 1 ]T (5) 

(x,y)£N(x,y) 

with a being the appropriate normalization coefficient and a neighbourhood N of (x,y). There are two 
important aspects to note here: 

• The resulting filter is by definition separable and therefore fast to compute. 


• We could also try O v j(Oh,i{J)) which would result in a new filter because Ohj and O v j are not 
commutative. An example can be seen in Figure [4j all neighbouring pixels of p get the same weight 
for both paths, except for q which gets a different weight. The filter Oh,i(O v j) does not detect a 
boundary between p and q whereas O v p(Ohj) uses a different path which has a boundary. 

In contrast to the 2D Gaussian range filter, this filter doesn’t directly compare I(x,y) and I(x,y) but 
rather uses an intermediate pixel I(x,y). Using the notation from above we have to = p = (x,p), 
t\ — (x, y) and £2 — <7 = (#, y)- The distance measure and constant that are necessary to define the path 
that we use here are d(u , v) = \u x — v x \ + | u y — v y \ and / = 1. 

To see how this separable filter differs from the 2D Gaussian range filter, we look at the computed weights 
in these two filters: 


,(p,g) ^ 


- 2 ^(/( p )-/(<7)) 2 


) 


W 1 ,(p,q) c 


(6) 

( 7 ) 


Analogous to [T| - [TV] this leads us to three different cases: 

• I(p) = I(ti) V I(q) = I(ti) w i,(p,q) = w 2 ,( P ,q) : If the intermediate pixel t\ equals p or g, then 
the resulting weight is the same as in the first approach. This corresponds to the cases [T| and [TT| 


• I(p) < I(ti) < I(q) V / (p) > I(ti) > I(q) => Wi^ p ^ > W 2 : (p, q )'- If the intermediate pixel t\ has 
an intensity value between p and q then although there may be a (large) difference between p and 
q (which the 2D Gaussian range filter would classify as an edge), this filter increases the resulting 
weight. That’s because although they differ, the intermediate pixel tells us they belong to the same 
region and that there’s a smooth transition from p to q. This corresponds to m 

• > max {/(p),/(g)}) V (/(U) < min {/(p),/(g)}) => w 1 ^ p ^ < w 2 ^p^ If the intermediate 
pixel t\ has a larger or smaller intensity value than both p and g, then there is an edge between p 
and q and the resulting weight is therefore smaller. This corresponds to m 

The path to, t \, t<i that we chose here leads to a reasonable description of a region for small neighbourhoods 
N(x,y). We therefore get a filter that is not only separable and therefore faster to compute than the 
2D Gaussian range filter, but that also uses an intermediate pixel to check if two pixels p and q he in 
the same region according to a path (see Figure [ 2 ] for an illustration of how the separable Gaussian filter 
works. Note how the 2D Gaussian range filter doesn’t recognize the edge between p and g3 whereas the 
separable Gaussian range filter correctly recognizes it.). 
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Figure 2: Schematic view of the 2D Gaus¬ 
sian range filter (left) and the separable Gaus¬ 
sian range filter (right). 



Figure 3: Schematic 
view of the SNN filter 
with a 5 x 5 neighbour¬ 
hood. 


Figure 4: Example 
showing that Ohj and 
O v j are not commuta¬ 
tive. 


3.3.4 Symmetric Nearest Neighbour Filter The third edge-aware filter we used is the symmetric 
nearest neighbour filter (SNN) first presented in [2]. It compares all pixels in a predefined neighbourhood 
of a pixel p and chooses of each pixel pair of opposite pixels the more similar neighbour to p. p is then 
replaced by the mean (SNNmean) or the median (SNNmedian) of the symmetric pixels more similar 
to the center. In the simplest case we consider a 3 x 3 neighbourhood, therefore it compares the eight 
neighbouring pixels of p and calculates the mean or median of four pixels. This leads to a very fast 
edge-aware filter. In our results section we considered only this simple 3x3 case. 


3.4 Advantages and Disadvantages 

The main advantage of this framework is that it is very flexible. According to the application field one 
can choose which filters are most suitable. Removing noise from an image may require a different guided 
filter than texture extraction or pre-processing for edge detection. 

Furthermore the SiR filter is very easy to implement. If you already have a smoothing and an edge- 
aware filter, then the implementation is achieved with only a few lines of code. The performance depends 
directly on the performance of both filters. Plugging an edge-aware filter into SiR will necessarily take 
more time than applying only the edge-aware filter alone, but it may improve the overall performance. 
There are however two drawbacks using this framework: First of all, as blurring may lead to less intense 
colors, this unwanted property transfers to the output image too, as the output image is a restoration 
based on the blurred image. While the resulting image consists of smooth regions and sharp edges, it 
may also lose some of its color intensity, especially if the whole algorithm is applied more than only 
one time. The second disadvantage concerns the restoration of small edges near large edges. Near large 
edges, which still exist after blurring, there exists color information to restore not only the large edge 
itself, but also surrounding structures with similar color, it is possible that small edges reappear near 
large edges after blurring them away. This means that near large structures, small structures or small 
artifacts can appear, because the edge-aware filter can partially restore them due to the available color 
information from the large edge. 

3.5 Comparison to the Rolling Guidance Filter 

RGF and SiR share the same main idea: remove small structures with an appropriate filter and use an 
edge-aware filter to handle strong edges. However they differ in how they implement this idea. The 
RGF fixes the input image, each iteration is filtering the same image. In every iteration it computes a 
new guidance image to better filter the input image in the next iteration, starting with a blurred input 
image in the first step. The SiR on the other hand does the reverse: it fixes the guidance image, using 
its optimal edge information, and instead modifies the image to be filtered in every iteration. It starts 
off with an initially blurred image and restores the strong edges according to the original input image. 
This differing setting of which image is fixed and which one changes iteratively leads to two different 
principles: calculate a reasonable guidance image to filter the input image once (RGF) versus use all the 
edge information available as guidance to restore important edges in a blurred image (SiR). In Section 
[4] you can see comparisons concerning the output of both algorithms. 

4 Applications and Results 

In this section we present results obtained with our framework during different experiments and show 
some possible applications. When calculating the Gaussian blur, the 2D Gaussian range filter or our 
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separable Gaussian range filter, we always used a 7 x 7 neighbourhood. For calculating the SNN we used 
a 3 x 3 neighbourhood. We applied our algorithm to each color channel separately in RGB color space. 
All results in this section where obtained on an Intel i7 processor with 2.9 GHz. Our C++ implementation 
uses no SIMD instructions and all implementations were running on a single thread, therefore there is still 
a lot of room for improvement. However, even this simple implementation results in a very fast filtering 
operation. The results are compared to the C++ implementation of the rolling guidance filter 
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using 

the joint bilateral filter which is freely available from their project homepage. It uses only one single 
thread and no SIMD instructions as well. However to speed up the bilateral filtering it implements the 
permutohedral lattice from 1 . 


4.1 Runtime evaluation 

Before we look at the actual output of our algorithm, we compare the runtimes of the rolling guidance 
filter (RGF) and different filter combinations using our model (SiR). In Table [l] we applied all listed 
filters to the image shown in Figure |5j For RGF there is only a total runtime entry as it doesn’t separate 
the smooth and restore steps. The parameters were chosen as follows: 

• RGF: a s = 3, a r = 25, 4 iterations 

• SiRSNN: SiR with iterated Box blur (5x5 neighbourhood, applied 2 times) and SNNmean. 

• SiR2DGauss: SiR with Gaussian blur (a = 5) and 2D Gaussian filter (a = 20) 

• SiRsep: SiR with Gaussian blur (a = 5) and separable Gaussian filter (a = 20) 

For SiR2DGauss and SiRsep we set the numbers of iterations to n m 5, for SiRSNN we set n = 9. 

As we can see, SiR yields far better runtime results than the rolling guidance filter, the box blur and 
SNN combination even results in lOfps real-time performance without any kind of optimization. Even 
our slowest method, SiR2DGauss, runs in only a third of the time of the RGF. 



RGF 

SiR2DGauss 

SiRsep 

SiRSNN 

Smooth step 

X 

0.055 

0.055 

0.035 

Iteratively Restore step 

X 

0.72s 

0.295 

0.07s 

Total 

2.4s 

0.775 

0.345 

0.15 


Table 1: Runtime comparison on the image seen in Figure 5 with dimensions 802x682 using no optimizations. 


4.2 Texture smoothing 

The smooth-and-iteratively-restore (SiR) framework can be used for texture smoothing or texture ex¬ 
traction. We tested it on various images and compared it to the rolling guidance filter (RGF), which 
is the only state-of-the-art scale-aware edge-preserving smoothing filter that we know of. All obtained 
results where very similar to those shown in Figure [5j In terms of smoothing regions with similar color, 
SiRsep achieves more homogeneous results than RGF and SiRSNN, which both produce bright stains. 
Furthermore the overall appearance of the edges with SiRsep seems better in terms of preserving the 
original shape of the edges. RGF produces smoother but also less detailed edges, while SiR tries generally 
to restore the edges as similar as possible to the original image. 

However, as mentioned in Section |3.4[ in the neighbourhood of edges SiRsep produces some artifacts 
because it partially restores small structures, while RGF doesn’t. Also, the resulting image of SiR suffers 
from less color intensity than the original one because of the smoothing step. 

SiRSNN produces noisier results than both SiRsep and RGF, however in terms of speed it is way ahead. 
If aiming for real time applications, SiRSNN may still satisfy your requirements while being extremely 
fast. 

If the output of RGF is more desirable, but the performance is too slow, one can combine our method 
with RGF: use SiR to compute a smoothed image with very good edge information fast and then use 
that image as initial guidance image in RGF and perform therefore only a few steps (ideally only one) 
with RGF. 
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Figure 5: Texture smoothing output of different filters 
SirGauss2D, SiRsep. Parameters are listed in Section 4.1 


, from left to right: original image, RGF, SiRSNN, 
The image was taken from 
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4.3 Pre-processing for edge detection 

Edge detection is a good example where small structures with strong edges can mislead the used edge 
detection algorithm. In such cases using a scale-aware edge-preserving filter can improve the result 
drastically by removing those small structures while retaining the important edges (for an example see 
Figure |6|. We tested our framework on the Berkeley Segmentation Dataset [2] together with the Sobel 
edge detector. We tested 3 different settings: edge detection without pre-processing, edge detection after 
using one of the three edge-aware filters and finally edge detection using those same edge-aware filters 
plugged in in our framework. As can be seen in Table [2] our algorithm does not only achieve significant 
better results than edge detection without pre-processing, it also achieves clearly better results than 
using one of the edge-aware filters only as pre-processing without the added scale-awareness provided by 
our method. For the edge-aware filters we used following parameters: 

• SNN: 3x3 neighbourhood 

• 2D Gaussian Range Filter: 7x7 neighbourhood, a = 8 

• Separable Gaussian Range Filter: 7x7 neighbourhood, a = 8 

When plugged into our framework we used the same parameters as above and additionally in each case 
we used a Gaussian blur in a 7 x 7 neighbourhood with a = 3 using 5 iterations. 



No filter 

SNN 

2D-Gauss 

sep. Gauss 

Without SiR 

0.49 

0.53 

0.51 

0.51 

With SiR 

X 

0.56 

0.56 

0.57 


Table 2: Benchmark results of Sobel edge detection using no pre-processing vs. edge-aware filter as pre-processing 
step vs. our method as pre-processing step on the Berkeley Segmentation Dataset The computed score represents 
the F-measure. 


5 Conclusion and further work 

In this paper we presented a novel framework for edge-preserving smoothing, which uses a smoothing 
filter for small structure removal and a guided edge-aware filter to restore strong edges. It can be used 
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Figure 6: Sample results of edge detection before and after applying SiRsep. Image taken from 4 . 


with a variety of different filters, of which a small subset was presented within this paper. The algorithm 
is very simple and therefore easy to implement. It consists of two steps: First smooth once using the 
smoothing filter, then iteratively apply the guided edge-aware filter. Depending on the performance of the 
two filters and the predefined number of iterations, the algorithm can be very fast. As our experiments 
show, with carefully chosen filters our framework is well suited for real-time applications. 

Further work on this topic includes experimenting with more guided edge-aware filters and applying the 
model to other applications. 
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